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The evidence for the accelerated expansion of the universe and the time-dependence of the fine- 
structure constant suggests the existence of at least one scalar field with a mass of order Ho. If such 
a field exists, then it is generally assumed that its coupling to matter must be tuned to unnaturally 
small values in order to satisfy the tests of the Equivalence Principle (EP). In this paper, we present 
an alternative explanation which allows scalar fields to evolve cosmologically while having couplings 
to matter of order unity. In our scenario, the mass of the fields depends on the local matter density: 
the interaction range is typically of order 1 mm on Earth (where the density is high) and of order 
10— 10 4 AU in the solar system (where the density is low). All current bounds from tests of General 
Relativity are satisfied. Nevertheless, we predict that near-future experiments that will test gravity 
in space will measure an effective Newton's constant different by order unity from that on Earth, 
as well as EP violations stronger than currently allowed by laboratory experiments. Such outcomes 
would constitute a smoking gun for our scenario. 



I. INTRODUCTION 

There is growing evidence in cosmology for the existence of nearly massless scalar fields in our Universe. On 
the one hand, a host of observations, from supernovae luminosity-distance measurements |l| to the cosmic microwave 
background anisotropy ,2] , suggests that 70% of the current energy budget consists of a dark energy fluid with negative 
pressure. While observations are consistent with a non-zero cosmological constant, the dark energy component is more 
generally modeled as quintessence: a scalar field rolling down a flat potential [3J, |^ . In order for the quintessence field 
to be evolving on cosmological time scales today, its mass must be of order Hq, the present Hubble parameter. 

On the other hand, recent measurements of absorption lines in quasar spectra suggest that the fine-structure 

n 

constant a has evolved by roughly one part in 10 5 over the redshift interval 0.2 < z < 3.7 yj. Time- variation of 



coupling constants are generally modeled with rolling scalar fields 



a. 



and the recent evidence for a time-varying a 



requires the mass of the corresponding scalar field to be of order Hq jj . 

In either case, the inferred scalar field is essentially massless on solar system scales, and therefore subject to tight 



constraints from tests of the Equivalence Principle (EP) |8|. The current bound on the Eotvos parameter, 77, which 
quantifies the deviation from the universality of free- fall, is r\ < 10 -13 , from the Eot-Wash experiment 

From a theoretical standpoint, massless scalar fields or moduli are abundant in string and supergravity theories. 
Indeed, generic compactifications of string theory result in a plethora of massless scalars in the low-energy, four- 
dimensional effective theory. However, these massless fields generally couple directly to matter with gravitational 
strength, and therefore lead to unacceptably large violations of the EP. Therefore, if the culprit for quintessence 
or time-varying a is one of the moduli of string theory, some mechanism must effectively suppress its EP-violating 
contributions. 

For instance, Damour and Polyakov |lo| (see also have proposed a dynamical mechanism to suppress the 

coupling constants (3i between the various matter fields and the dilaton of string theory. Alternatively, the suppression 



could be the result of approximate global symmetries . 

In a recent paper , we presented a scenario in which scalar fields can evolve cosmologically while having couplings 
to matter of order unity, i.e., /3, ~ This is because the scalar fields acquire a mass whose magnitude depends 

on the local matter density. In a region of high density, such as on Earth, the mass of the fields is large, and thus the 
resulting violations of the EP are exponentially suppressed. In the solar system, where the density is much lower, the 
moduli are essentially free, with a Compton wavelength that can be much larger than the size of the solar system. 
Finally, on cosmological scales, where the density is very low, the mass can be of the order of the present Hubble 
parameter, thereby making the fields potential candidates for causing the acceleration of the universe or the time- 
evolution of the fine-structure constant. While the idea of density-dependent mass terms is not new 
the novelty of our work lies in that the scalar fields can couple directly to baryons with gravitational strength. 

In our scenario, scalar fields that have cosmological effects, such as quintessence, do not result in large violations of 
the EP in the laboratory because we happen to live in a very dense environment. Thus, the main constraint on our 
model is that the mass of the field be sufficiently large on Earth to evade EP and fifth force constraints ^1 ■ 

The generation of a density-dependent mass for a given modulus <f> results from the interplay of two source terms in 
its equation of motion. The first term arises from self-interactions, described by a monotonically-decreasing potential 
V(4>) which is of the runaway form (see Fig. In particular, we underscore the fact that the potential need not 
have a minimum; rather, it must be monotonic. The second term arises from the conformal coupling to matter fields, 
of the form e^' Mpi . The coupling constants Pi need not be small, however, and values of order unity or greater 
are allowed. Although these two contributions are both monotonic functions of <j), their combined effect is that of 
an effective potential which does display a minimum (see Fig. [5J. Furthermore, since this effective potential depends 
explicitly on the local matter density p, both the field value at the minimum and the mass of small fluctuations depend 
on p as well, with the latter being an increasing function of the density. 

Although the scalar fields are quite massive on Earth, their behavior is strikingly different in the solar system where 
the local matter density is much smaller. Thus our model makes a crucial prediction for near-future experiments 
that will test gravity in space. For example, consider the SEE Project which among other things will measure 
Newton's constant to an unprecedented accuracy. Our scenario generically predicts that the SEE experiment should 
observe corrections of order unity to Newton's constant compared to its measured value on Earth, due to fifth-force 
contributions which are important in space but exponentially suppressed on Earth. 

Moreover, three satellite experiments to be launched in the near future, STEP 19], Galileo Galilei (GG) 20j and 
MICROSCOPE 0, will test the universality of free-fall in orbit with expected accuracy of 1CT 18 , 1CT 17 and 1CT 15 , 
respectively. We predict that these experiments should observe a strong EP- violating signal. In fact, for a wide range 
of parameters, our model predicts that the signal will be larger than the ground-based Eot-Wash bound of 10 -13 . 

If SEE does measure an effective Newton's constant different from that on Earth, or if STEP observes an EP- 
violating signal larger than thought permitted by the Eot-Wash experiment, this will strongly indicate that a mech- 
anism of the form proposed here is realized in Nature. For otherwise it would be hard to explain the discrepancies 
between measurements in the laboratory and those in orbit. These new and surprising outcomes are a direct conse- 
quence of the fact that scalar fields in our model have drastically different behavior in regions of high density than in 
regions of low density. 

We refer to <fi as a "chameleon" field, since its physical properties, such as its mass, depend sensitively on the envi- 
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ronment. Moreover, in regions of high density, the chameleon "blends" with its environment and becomes essentially 
invisible to searches for EP violation and fifth force. 

Even though we predict significant violations of the EP in space, all existing constraints from planetary orbits Q], 
such as those from lunar laser ranging are easily satisfied in our model. This is because of the fact that the 
chameleon-mediated force between two large objects, such as the Earth and the Sun, is much weaker than one would 
naively expect. To see this, we use calculus and break up the Earth into a collection of infinitesimal volume elements. 
Consider one such volume element located well- within the Earth. Since the mass of the chameleon is very large inside 
the Earth, the 0-flux from this volume clement is exponentially suppressed and therefore contributes negligibly to the 
0-ficld outside the Earth. This is true for all volume elements within the Earth, except for those located in a thin shell 
near the surface. Infinitesimal elements within this shell are so close to the surface that they do not suffer from the 
bulk exponential suppression. Thus, the exterior field is generated almost entirely by this thin shell, whereas the bulk 
of the Earth contributes negligibly. A similar argument applies to the Sun. Consequently, the chameleon-mediated 
force between the Earth and the Sun is suppressed by this thin-shell effect, which thereby ensures that solar system 
tests of gravity are satisfied. 

However, note that this only applies to large objects, such as planets. Sufficiently small objects do not suffer from 
thin-shell suppression, and thus their entire mass contributes to the exterior field. In particular, a small satellite in 
orbit, such as SEE, may not exhibit a thin-shell effect. This is why the orbits of the planets are essentially unaffected 
by the 0-force, whereas the fifth force between two test particles in the SEE capsule is significant. 

Since <fi couples directly to matter fields, all mass scales and coupling constants of the Standard Model depend on 
space and time. Once again due to the thin-shell mechanism described above, spatial variations of constants are suffi- 
ciently small in our model to satisfy current experimental bounds, for instance from the Vessot-Levine experiment [2^- 
Time variation of coupling constants are also not a problem since, during most of the history of the universe, the 
various couplings actually change by very little. Thus the bounds from big bang nucleosynthesis, for instance, are 
easily satisfied. This will be described in more detail in a separate paper dealing with the cosmological evolution in 
our model |2J] . 

In Sec.[nl we describe the ingredients of the scenario, focusing on a single modulus (j> for simplicity. We show how 
the dynamics of <f> are governed by an effective potential that depends on the local matter density. In Sec. IIIII we 
derive approximate solutions for <p for a compact object such as the Sun, for instance, and describe the thin-shell 
mechanism mentioned earlier. In Sec. IIVI we specialize the solution for <f> to the case of the Earth, and apply the 
results in Sec.0to derive constraints on the parameters of the theory based on laboratory tests of the EP and searches 
for a fifth force. We then show in Sec. lVII that. for a potential of power-law form, V((f>) = M 4+n (j>~ n , these constraints 
translate into the requirement that the energy scale M be less than an inverse millimeter or so. Curiously, this is 
also the scale associated with the cosmological constant today. In Sec. IVIll we argue that our model easily satisfies 
constraints from solar system tests of GR. It is showed (Sec. IVIIIf) that the same holds true for bounds from spatial 
and time variation of coupling constants. We then predict (Sec. IIXII that near-future experiments that aim at testing 
the EP and measuring a fifth force should observe a large signal, perhaps stronger than previously thought possible. 
Finally, we conclude and summarize our results in Sec. Ixl 
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FIG. 1: Example of a runaway potential. 

II. THE INGREDIENTS OF THE MODEL 

Focusing on a single scalar field <f> for simplicity, the action governing the dynamics of our model is given by 

S= [ d'xV—gl^n-Udtf -V(A - [ d 4 z£ m (^\ 5 «), (1) 
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where Mpi = is the reduced Planck mass, g is the determinant of the metric g M „, TZ is the Ricci scalar and 

(i) 

ipm are matter fields. The scalar field <p interacts directly with matter particles through a conformal coupling of the 
form e^^ Mpi . Explicitly, each matter field ipm couples to a metric gfy which is related to the Einstein- frame metric 
9fiv by the rescaling 

g(i) =e ^/M Plg ^ : ( 2 ) 

where /3j are dimensionless constants 25] . Moreover, for simplicity, we assume that the different ipm 's do not interact 
with each other. Note that Eq. JQ) is of the general form of low-energy effective actions from string theory and 
supergravity, where V(<p) arises from non-perturbative effects. 

The potential V(<p) is assumed to be of the runaway form. That is, it is monotonically decreasing and satisfies 

lim V = , lim ^ = , lim = . . . , (3) 

as well as 

lim V = oo , lim — — = oo , lim — = oo . . . , (4) 
where = dV/dcp, etc. See Fig. ^ For instance, a fiducial example is the inverse power-law potential 

V{<j)) = M 4+n <p~ n , (5) 

where M has units of mass and n is a positive constant. The above conditions on the asymptotics of V are generally 
satisfied by potentials arising from non-perturbative effects in string theory [lj, lift lid \2(\ . Note that it is also of the 
desired form for quintessence models of the universe [2^ . 
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The equation of motion for <f> derived from the above action is then 

where T$ = (2/y/— gM)5C m /Sg^ is the stress-energy tensor for the ith form of matter. For the purpose of this 
paper, it will suffice to approximate the geometry as Minkowski space, that is, g^ v ~ rj^. This is valid provided that 
the Newtonian potential is small everywhere, and that the backreaction due to the energy density in <fi is also small. 
This latter assumption will be justified when we analyze post-Newtonian corrections in Sec. IVII Bl 

For non-relativistic matter, one has g^T^l! w ~Pi, where pi is the energy density. However, we shall find it 
convenient to express our equations not in terms of pi, but rather in terms of an energy density pi = p i e^ i ^' Ml>l 
which is conserved in Einstein frame. In other words, pi is defined so that it is independent of cj>. Equation © thus 
reduces to 

v 2 ^^ + £ifW I0/Mpi - (7) 

. Mpl 

I 

From the right-hand side of Eq. J7J), we see that the dynamics of (p are not solely governed by V ((/)), but rather by 
an effective potential 

V eff {cp)^V{cp)+Y.P^ /Mpl ( 8 ) 

i 

which depends explicitly on the matter density pi. In particular, although V((f>) is monotonic, V e ff does exhibit a 
minimum provided that /3j > 0. This is illustrated in Fig. [21 for the case of a single component p with coupling (3. 
(One could equivalently consider the case V^, > and (3i < 0.) Unfortunately, known examples in string theory have 
and Pi occurring with the same sign. If 4> is the modulus of an extra dimension, for instance, then one expects 
that V — > and — > in the decompactification limit in this limit, the various masses of particles also tend 
to zero. 

We will denote by (j)min the value assumed by (j) at the minimum, that is, 

V^ mm ) + V J^Lpje/fc = . (9) 

*-r Mpi 

Meanwhile, the mass of small fluctuations about 4>min is obtained as usual by evaluating the second derivative of the 
potential at 4> m in'- 

mlin = VMM + E -§ rPl e^-' M r< . (10) 

M Pl 



Equations © and (|10f) respectively imply that the local value of the field, 4> m im and its mass, m m i n , both depend 
on the local matter density. Since is negative and monotonically increasing, while is positive and decreasing, 
it also follows that larger values of pi correspond to smaller 4> m in and larger m m in- This is illustrated in Fig. 
The denser the environment, the more massive is the chameleon. We will later see that it is possible for m m i n to be 
sufficiently large on Earth to evade current constraints on EP violations and fifth force, while being sufficiently small 
on cosmological scales for (j> to have interesting cosmological effects. 

The upshot of our model, from a theoretical standpoint, is that the potential V(<p) need not have a minimum, nor 
does the coupling constant j3 need be tuned to less than 1CP 4 to satisfy EP constraints ^^|. Quite the contrary, V{4>) 
is assumed monotonic, while (3 can be of order unity. 
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FIG. 2: The chameleon effective potential V e ff (solid curve) is the sum of two contributions: one from the actual potential 
V(<f>) (dashed curve), and the other from its coupling to the matter density p (dotted curve). 
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Large p Small p 

FIG. 3: Chameleon effective potential for large and small p, respectively. This illustrates that, as p decreases, the minimum 
shifts to larger values of <j) an d the mass of small fluctuations decreases. (Line styles are the same as in Fig. |5|) 

III. PROFILE FOR A COMPACT OBJECT 

In order to study the observable consequences of our model, in particular with regards to EP violations and fifth 
force mediation, we must first understand the profile that (f> acquires on Earth and in the solar system. Therefore, 
in this Section, we derive an approximate solution for <j) in the case where the source is a compact object, which we 
idealize as being perfectly spherical and having homogeneous density. 

Thus consider a static, spherically-symmetric body of radius R Cl homogeneous density p c and total mass M c — 
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4ttR^Pc/3. In the case of the Earth, the latter can be approximated by the characteristic terrestrial density: p® w 
10 g/cm 3 . We treat the object as isolated, in the sense that the effect of surrounding bodies is neglected. It is 
not in vacuum, however, but is instead immersed in a background of homogeneous density p^ . In the case of solar 
system objects, this models the fact that our local neighborhood of the galaxy is not empty, but rather filled with 
an approximately homogeneous component of baryonic gas and dark matter with density poo = pa ~ 10~ g/cm 3 . 
In the case of a baseball in the Earth's atmosphere, p^ denotes the surrounding atmospheric density: = p a tm ~ 
10" 3 g/cm 3 . 
With these assumptions, Eq. (JJJ reduces to 

3 <"> 

where 

*)-{» f " r< V (12) 
ypoo for r > R c 

Note that we temporarily focus on the case where all /3j's assume the same value (3. This is done for simplicity 
only, and the following analysis remains qualitatively unchanged when these are taken to be different. Moreover, this 
assumption will be dropped when we derive the resulting violations of the EP in Sec.0 In other words, in the end 
we are not assuming that the theory is Brans-Dicke j^l ■ 

Throughout the analysis, we denote by 4> c and 4>oo the field value which minimizes V e ff for r < R c and r > R Cl 
respectively. That is, from Eqs. © and l(T2|> . we have 

Mpi 

K^oo) + -^-Pooe^' Mpl = . (13) 
Mpi 

Similarly, we denote by m c and the mass of small fluctuations about <j) c and </>oo, respectively. That is, m c (moo) 
is the mass of the chameleon field inside (outside) the object. 

Equation is a second order differential equation and as such requires two boundary conditions. Since the 
solution must be non-singular at the origin, we require d<j>/dr — at r = 0, as usual. Moreover, since p = at 
infinity, it is natural to impose cj> — > 4>oo as r — > oo. This latter condition is physically sensible as it implies dtf>/dr — ► 
as r — > oo, and thus that the </)-force between the compact object and a test particle tends to zero as their separation 
becomes infinite. To summarize, Eq. Ijlljl is subject to the following two boundary conditions: 

^ = at r = 0; 
dr 

4> ^ 4>oo as v — > oo . (14) 



A. Qualitative description of the solution 



Before solving this problem explicitly it is useful to give a heuristic derivation of the solution. Far outside the 
object, r 3> R c , we know that the chameleon tends to </>oo, as required by the second boundary condition. There are 
then two types of solution, depending on whether the object is large {e.g., the Earth) or small (e.g., a baseball). The 
distinction between large and small will be made precise below. 
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FIG. 4: For large objects, the 0-field a distance r > R c from the center is to a good approximation entirely determined by the 
contribution from infinitesimal volume elements dV (dark rectangle) lying within a thin shell of thickness AR C (shaded region). 
This thin-shell effect suppresses the resulting chameleon force. 



Small objects do not generate large variations in <f). Thus, their solution can be thought of as a perturbation on the 
background solution = 0oo. Hence, one has 4> « everywhere in this case, including the interior of the object. 

Large objects, on the other hand, are strongly perturbing. Within the object, r < R Cl one finds that the chameleon 
nearly minimizes the effective potential, and thus (f> ~ <j) c . Hence, the solution essentially extrapolates between <f> — <f> c 
within the core and <fi = </>oo far outside. 

To be more precise, let us describe the exterior solution (r > R c ) for large bodies, assuming m^Rc -C 1 for simplicity. 
For this purpose, it is convenient to break up the object into infinitesimal volume elements dV and consider their 
individual contribution to the ^-field. Well-within the object, one has <p ~ <p c , and the mass of the chameleon is large, 
m c 3> rritx,. Thus, the contribution from a volume element dV within the core is proportional to exp(— m c r) and is 
therefore exponentially suppressed. In other words, it contributes negligibly to the </>-field outside. This holds for all 
infinitesimal volume elements within the object, except for those lying within a thin shell of thickness AR C near the 
surface |3Cl|. See Fig. 0] Thus, the exterior solution is obtained by summing over all elements within this shell: 



(16) 



^ttMpi J \ R. 

where r is the distance from the center of the object. We will find in Sec. IIIIBI that 

R c ~ Qf3Mp&c ' 

where $ c = Mc/S'kM^Rc is the Newtonian potential at the surface of the object. The derivation of Eq. (|15fl implicitly 
assumed that the shell was thin, that is, 

AR C 

R c 

We shall henceforth refer to this as the thin-shell condition 



<1. (17) 
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Small objects, in the sense that AR C /R C > 1, do not have a thin shell. Rather, their entire volume contributes to 
the 0-field outside, and thus the exterior solution is 



\4ttMpi 

which is recognized as the Yukawa profile for a scalar field of mass irioo . Note that Eq. H15|) and H18|) differ only by a 
thin-shell suppression factor of AR C / R c . To summarize, the exterior solution for a compact object is given by 

( P \f3AR A M c e-^ _ if ^£ <<1) (19) 



\AttMpiJ V, Rc J r R c 
with AR C /R C defined in Eq. fity . 

The ratio ((poo — (p c )/Mpi<f> c , which appears in Eq. (|f 6() and determines whether or not an object has a thin shell, 
can be interpreted physically as follows. Given a background profile (p(r), it is straightforward to show from the 
action that the resulting chameleon-force on a test particle of mass M and coupling f3 is given by 

F = --|-MW. (20) 

Mpi 

It follows that (p should be thought of as a potential for this "fifth" force. Thus, ((poo — (p c )/Mpi& c is the ratio of the 
difference in <fi potential to the Newtonian potential, and effectively quantifies how perturbing the object is for the <p 
field. 



B. Derivation 

To get an intuition for the boundary value problem at hand, it is useful to think of r as a time coordinate and 
<p as the position of a particle, and treat Eq. ({11(1 as a dynamical problem in classical mechanics. This is akin to 
the familiar trick performed in bubble nucleation calculations 31]. In this language, the particle moves along the 
inverted potential —V e fj, and the second term on the left-hand side of Eq. (|llf> . proportional to 1/r, is recognized 
as a damping term. An important difference here is that — V e ff is "time" -dependent since p depends on r. More 
precisely, the effective potential undergoes a jump at time r = R c , as illustrated in Fig. El 

The particle begins at rest, since d<p/dr = at r = 0, from some initial value which we denote by (pi'- 

fa ee (P(r = 0) . (21) 

For small r, the friction term is large, and thus the particle is essentially frozen at <p = (pi. It remains stuck there 
until the damping term, proportional to 1/r, is sufficiently small to allow the driving term, dV e ff /d(p, to be effective. 
In other words, the amount of "time" the particle remains stuck near (p = fa depends on the slope of the potential, 
dV e ff/d(p, at (p = fa. Once friction is negligible, the particle begins to roll down the potential. See Fig. EK) . 

It rolls down until, at some later time r = R c , the potential suddenly changes shape as p(r) undergoes a jump from 
Pc to poo (see Eq. I|12|l ). But (p and d(p/dr are of course continuous at the jump, and the particle keeps rolling, this 
time climbing up the inverted potential. See Fig. \Ejp) . If the initial position fa is carefully chosen, the particle will 
barely reach (poo in the limit r — > oo, as desired. (It is easy to prove that such fa always exists.) Thus the problem is 
reduced to determining the initial value fa. 
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p inf exp(p>/ Mpj) 



p exp(p>/ M ) 
c 





-v (4>) 

eff v T/ 



-V ((j)) 

eff v T/ 



a) r < R c 



b) r > R c 



FIG. 5: The inverted potential — V e // for a compact object of radius i? c is discontinuous at r — R c since the matter density 
equals: a) p c for r < R c ; b) for r > R c . The dots represent the position of the particle at some value of r. 



In the end fa will depend on the physical properties of the compact object, such as its density p c and its radius 
R c , as well as on the various parameters of the theory, such as (3 and the shape of the potential. Rather than first 
choosing a set of values for p c , R c , etc. and then solving for fa, we will instead choose a range of fa and determine 
the corresponding region in the (p c , R c , . . .) parameter space. More precisely, we shall consider the two regimes 
(fa — (j> c ) <C fa and fa > fa, and we will show that these correspond respectively to AR C /R C <C 1 and AR C /R C > 1. 
Anticipating this result, we refer to these two regimes as thin-shell and thick-shell, respectively. 

Thin-shell regime: (fa — fa) <C fa. This corresponds to cj) being released from a point very close to fa. Since <f> c 
is a local extremum of the effective potential, the driving term dV e ff/d(f> is negligible initially, and the dynamics are 
strongly dominated by friction. Consequently, the field remains frozen at its initial value fa w fa for a long time, 
until the friction force is sufficiently small to allow the particle to roll. We shall denote by R ro u the "moment" at 
which this occurs. Hence, we have 



(r) w fa for < r < R roll 



(22) 



When r ~ R ro ii, the field is still near fa but has now begun to roll. Since Mpi \V,p\ <C ^pe^^l Mpl as soon as cf> is 
displaced significantly from fa, as illustrated in Fig.|3 we may approximate Eq. 1)11(1 in the regime R ro u < r < R c 
by 

fl + l^l- JL pc (23) 
dr 2 r dr Mpi ' 

where we have also assumed j3<f>/Mpi -C 1. The solution to Eq. (|23(l with boundary conditions <fi = 4>c and dip/dr = 
at r = Rroii is 



— + 



R 



roll 



PPcRloh 



for R ro a <r < R c 



(24) 



Ml pi V 2 r J 2Mpi 
The full solution for < r < R c is thus approximated by Eqs. (|22|) and l|24|) . The approximation of separating the 
solution into the two regions < r < R ro ii and R ro ii < r < R c only makes sense, however, if R c — R ro ii "C R c - For 



11 



otherwise there is no clear separation between the two regions, and one needs a solution valid over the entire range 
< r < R c . 

At r = R c , the energy density undergoes a jump from p c to p^,, as described by Eq. (|12[1. For r > R c , the effective 
potential is shown in Fig.[5]b), and the particle is climbing up the hill. Since the speed of the particle is initially large 
compared to the curvature of the potential, Eq. (jl 1|> can be approximated by 

^ + -^«0 (25) 
dr 2 r dr 



whose solution, satisfying 4> — * 4>oo as r — > oo, is given by 

(7 e -"ioo (r-Ra) 



+ <f>oo , (26) 



where C is a constant and where we have used the fact that the potential is approximately quadratic near (j) — . 

The two unknowns, R ro ii and C, are then determined by matching <f> and dtfi/dr at r = R c using Eqs. (|24() and 12611 . 
With the approximation that R c — R ro u <C i? c , it is straightforward to show that the exterior solution is 



« 1 , (28) 



An M pi J V R, 

with 

AR C _ — 4> c ^ R c — Rroll 

~rT = 6/3M Pi $ c ~ R c 
where we have substituted the Newtonian potential $ c . 

Thick-shell regime: <fri > <fr c . In this case, the field is initially sufficiently displaced from <p c that it begins to roll 
almost as soon as it is released at r — 0. Hence there is no friction-dominated regime in this case, and the interior 
solution for <j> is most easily obtained by taking the R ro ii limit of Eq. (|24Jl and replacing <f> c by 4>i . Matching to 



the exterior solution as before, we obtain 



QMpi 

and 



0(r) = +(f>i for < r < R c (29) 



/ 3 \ M e _m °°( r_ - Rc ) 
(r) w - --— — + (^oo for r>R c . (30) 



Moreover, equating these two equations at r — R Cl we find <f>i = — 3/3$ c /Mp;. In particular, since fa > <p c and 
using the definition AR C /R C — (cfroo — c )/6/3A/p;$ c , this implies 

AR C 



R, 



>1. (31) 



We conclude this Section with a word on how the above solutions, which assumed homogeneous p, can be generalized 
to the more realistic case of spatially- varying matter density. In most cases of interest, such as the interior of the 
Earth for instance, we will find that the matter density varies on scales much larger than the Compton wavelength 
ui~ of the chameleon field in that region. More precisely, it is generally the case that |V \ogp(x)\ <C m within dense 
objects. If so, one can make an adiabatic approximation which consists of treating p(x) as a constant in the equations 
of motion. In other words, in this case one may simply substitute p(x) in the expressions above. 
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IV. PROFILE FOR THE EARTH 

Since the most stringent constraints on possible violations of the EP derive from experiments performed on Earth, it 
is important to discuss in some detail the profile for <p inside and in the vicinity of the Earth. Admittedly, the model 
for our planet described below is rather crude, but is sufficiently accurate to derive order-of-magnitude estimates 
of resulting violations of the EP. More realistic descriptions can be obtained for instance by using the adiabatic 
approximation discussed at the end of Sec. II I II or through numerical analysis. 

The Earth is modeled as a solid sphere of radius R® — 6 • 10 8 cm and homogeneous density = 10 g/cm 3 . 
Surrounding it is an atmosphere which we approximate as a layer 10 km in radius with homogeneous density p a tm — 
10 -3 g/cm 3 . Moreover, we treat our planet as an isolated body, neglecting the effect of surrounding compact objects 
such as the Sun and the Moon. Furthermore, far away from the Earth, the matter density is approximated by the 
density of homogeneous gas and dark matter in our local neighborhood of the galaxy: pc = 10 -24 g/cm 3 . 

The set-up is thus almost identical to that of Eq. Qlljl. except that the matter density now has three phases instead 
of two: 



p{r) 



p e for < r < i? ffi 

Patm for R® <r < R atm , (32) 
Pg for r > Ratm 



where R a tm = -R®+ 10 km. We henceforth denote by 0®, 4> a tm and 4>g the field value which minimizes the effective 
potential for the respective densities. Similarly, m®, m a tm and tug are the respective masses. 

Following the discussion in Sec. II I II the solution depends on whether or not the Earth and its atmosphere have a 
thin shell. As we will prove below, it is necessary that the atmosphere has a thin shell, for otherwise unacceptably 
large violations of the EP will ensue. In this case, one has cj> « 4> a tm in the bulk of the atmosphere. Moreover, since 
the Earth is much denser than the atmosphere, it follows that the Earth itself has a thin shell, in which case (f> ~ 0® 
inside the Earth. 

From Eqs. Ijl6(l and (|17|l . the thin-shell condition for the atmosphere reads 

ARatm _ $G — (j>atm ^ 
Ratm ~ WMp^atm ' 1 ' 

where $ atm = /9 atm i? 2 tm /6M|, ; . We can refine this bound by noting that, in order for the atmosphere to have a thin 
shell, clearly the thickness of the shell must be less than the thickness of the atmosphere itself, which is ~ 10~ 3 R a tm- 
Hence this requires AR at m/Ratm 10~ 3 . Using the fact that p atm ~ 10 _4 /O®, and thus $„ tra « 10 _4 $®, we can 
write this as 

A-B e _ 4>G ~ <t>atm »_ 7 .-.v 

R @ - 6/3M Pi $ e ' 1 ' 

This condition, which ensures that the atmosphere has a thin shell, will play a crucial role in the analysis of tests 
of gravity in the following sections. The exterior solution, r > Ratm, is then given by Eq. Q27J) with moo = ttig and 

0OO = <pG- 
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To summarize, the solution for <fi for the Earth and its atmosphere is well-approximated by 

4>{r) » 4>(s for < r < R B ; 

4>(r) » 0atm for i?© < r < i? atm ; 

0(r) w - - af — — h G for r > i? a tm , (36) 



V47rAf Pi J \ Rq 

where Ai? ffi /i? e is defined in Eq. (|3IJ> . 

It remains to show that tests of the EP require the atmosphere to have a thin shell. The proof proceeds by 
contradiction. Suppose that condition l|34|l is violated and instead we have 

^®>1(T 7 . (37) 

R<s 

We can therefore ignore the atmosphere altogether, and thus the </>-profile in the laboratory is given by Eq. (|35|l : 

where we have neglected the exponential factor since mcR® -C 1, as we will see in Sec. lVIl From Eq. Q20|l. this profile 
results in a fifth force on a test particle of mass M and coupling /3j of magnitude 

Supposing that the fy's are all of order (5 but assume different values for different matter species, then the resulting 
difference in relative free-fall acceleration for two bodies of different composition will be 

^2^-^-10-/^, (40) 

Q, 

as used in the Eot-Wash experiment 9j. For (3 of order unity, we see from Eq. (|37|) that this violates the bound 
r) < 10~ 13 . It follows that the atmosphere must be have a thin shell. 

V. SEARCHES FOR EP VIOLATION AND FIFTH FORCE ON EARTH 

The tightest constraints on our model derive from laboratory tests of the EP and searches for a fifth force HQ 
Since these experiments are usually done in vacuum, we first need to derive an approximate solution for the chameleon 
inside a vacuum chamber. For simplicity, we model the chamber as a perfectly empty, spherical cavity of radius R vac - 
In the absence of any other parts within the chamber, and ignoring the effect of the walls, the equation for <j) is given 
by Eq. CEH): 

d 2 2d<j> 13 

-d^ + -rTr =V ^ + M7 l p[r) > (41) 
where we have assumed (5<j)/Mpi < 1 as we did throughout Sec. IIIII and where 

for r < R vac 

(42) 

for r > R vac 



where 77 is the Eotvos parameter, and where the numerical coefficient is appropriate for Cu and Be test masses 
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The boundary conditions are the same as in Sec. IIIII dtfi/dr — at r — and <j» — > <j) a tm as r -> oo. 

The solution within the vacuum chamber is analogous to the solution for a compact object with thin-shell (see 
Sec. In both cases, due to the large density contrast between the object or the vacuum cavity and their 

environment, <p must start at r = from a point where it can remain almost frozen for the entire volume. That is, 
for < r < R c in the case of the overdense object and < r < R vac for the vacuum cavity. For a compact object, 
this freezing point lies naturally near the local extremum <fi = <fi c of the effective potential. For the vacuum chamber, 
the effective potential has no extremum for r < R vac (since p = there), and thus the only way <p can remain still is 
by starting from a point where it will be friction-dominated for almost the entire range < r < R vac - In other words 
the curvature of the potential at that point must be of order -R^" 2 c , that is, this freezing point corresponds to a value 
4> = 4>vac where the mass of small fluctuations, m vac , is equal to R^ac- 

While the precise solution to lj4*T|) depends of course on the details of the potential, numerical analysis confirms the 
above qualitative discussion: 

• The chameleon assumes the value <j) ~ <t>vac within the vacuum chamber, where <fi vac satisfies 

m vac = yV^f^w) = R~ac ■ (43) 

That is, 4>vac is the field value about which the Compton wavelength of small fluctuations equals R vac , the 
radius of the chamber. 

• Throughout the chamber, <f> varies slowly, with \dcf>/dr\ < (j) vac / 1 R va c- 

• Outside the chamber the solution tends to <fi a tm within a distance of m~^ m from the walls. 
These generic properties are all we need for our analysis. 



A. Fifth Force Searches 



The potential energy associated with fifth force interactions is generally parameterized by a Yukawa potential: 

Mi Mo e~' r/A 

V(r)=-a 1 „ , (44) 

where M\ and M% are the masses of two test bodies, r is their separation, a is the strength of the interaction (with 
a = 1 for gravitational strength), and A is the range. Null fifth- force searches therefore constrain regions in the (A, a) 
parameter space (see Fig. 2.13 of 

As discussed above, the range A of (/(-mediated interactions inside a vacuum chamber is of the order of the size of 
the chamber. That is, A ~ Rvac- For A « 10 cm — 1 m, the tightest bound on the coupling constant a from laboratory 
experiments is from Hoskins et al. [s^ : 

a < 1CT 3 . (45) 

Now consider two identical test bodies of uniform density p c , radius R c and total mass M c . If these have no thin 
shell, then they each generate a field profile given by Eq. (|18fl with (f)^ = 4> vac and moo = Rvac- 

^"-(Sdwfr 1 **- (46) 
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Dropping the irrelevant constant, the resulting potential energy is 

Comparison with Eq. I|44|l shows that the coupling strength is a = 2/3 2 in this case, which clearly violates the bound 
in Eq. IpES) for 8 ~ 0(1). 

Hence it must be that the test masses have a thin shell, that is, 

AR C _ <p vac - <f) c 

In this case, their field profile is given by Eq. Ijl5|) : 

(3 \ /3Ai? c \ M c e- r / R — 



AttMpi J \ R 



(49) 



and the corresponding potential energy is 



pi 



V R. 

Once again comparing with Eq. Q44JI. we find that the bound in Eq. I|45f) translates into 

2/32 (^r) 2 ~ 10-3 • (51) 

Note that, (3 > O(l), this constraint implies that the thin-shell condition in Eq. I|48|) is satisfied. 

To make the condition l|51|l more explicit, note that a typical test body used in Hoskins et al. had mass M c m 40 g 
and radius R c ~ 1 cm, corresponding to $ c « 3 • 10~ . Substituting in Eq. (|48|l and assuming <fi vac ^ <fi c , we obtain 
the constraint 

<t>vac < 10- 2S Mpi , (52) 
which ensures that the current bounds from laboratory searches of a fifth force are satisfied. 

B. Tests of the EP 

Turning our attention to the magnitude of EP violations inside our vacuum cavity, we recall that the solution for 
<i> inside the chamber satisfies 



d(i> 
dr 



(53) 



Rvac 

From Eq. H2U|I . this yields an extra component to the acceleration of magnitude (J3 /Mpi)cf) vac / R vac . For Cu and Be 
test masses, as used in the Eot-Wash experiment |9j, this yields an Eotvos parameter of 

r,~4nW-*p M rf® ^L. (54) 
Substituting R vac = 10 cm and /3 ~ 0(1), the Eot-Wash constraint of ij < 10 -13 translates into 

4>vac < 10~ 26 M Pl , (55) 
which is a much weaker constraint than Eq. il~>21 and thus shall henceforth be ignored. 
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VI. RESULTING CONSTRAINTS ON MODEL PARAMETERS 

In this Section we summarize the constraints derived in the previous sections and apply them to a general power-law 
potential 

V{<j)) = M 4+n 0~ n , (56) 

where M has units of mass and n is a positive constant. As mentioned earlier, potentials of this form have the desired 
features for quintessence models of the universe j^J. We will find that the energy scale M is generally constrained 
to be of the order of (1 mm) -1 . We then discuss the resulting bounds on the interaction range in the atmosphere, in 
the solar system and on cosmological scales. 

Broad considerations of EP violation lead us to conclude in Sec. HVI that both the Earth and its atmosphere must 
have a thin shell. This resulted in the constraint (see Eq. 



Ai? e _ 4> G 



<1(T 7 . (57) 



Recall that, by definition, 4>g is the value of cj> which minimizes the effective potential with p = pa, i.e., V^(<Pg) + 
(3pG^^ G l Mpl /Mpi = 0. Assuming j3<pG/Mpi <C 1 and substituting the power-law potential of Eq. l|56|l gives 

iM*+"M Pl \ 1/{n+1) 



PPG 

With pa = 1CP 24 g/cm 3 and <I>® = 10~ 9 , it is then straightforward to show that Eq. I|57|l can be rewritten as a bound 



/ 6 n+1 \ "+ 4 n+2 lBn-7 

M < ft—* ■ IO^ht . (l mm)" 1 . (59) 



on M: 



n 

Then, in Sec. |VJ we studied laboratory tests of the EP and fifth force, including the fact that these are performed 
in vacuum, and derived the condition 

4>vac < 1(T 28 M Pl , (60) 

where 4> vac is the field value about which the Compton wavelength of small fluctuations is of the order of the size of 
the vacuum chamber, R vac . Applying Eq. (|43|) to the power-law potential yields 

R~l = 1(T 68 M Pl = n(n + l)M 4 +>-j™+ 2 ) , (61) 

where we have assumed R vac = 1 m = 10 34 Mp, 1 for concreteness. It is easy to sec that smaller values of R vac result 
in weaker constraints, and this is therefore a conservative choice. It is then straightforward to show that Eq. 152|) 
reduces 

M < [n{n + i)]-V(4+") . io 3 ™/( 4 +") . (1 mm)- 1 . (62) 

We thus see that, for n and (3 of order unity, Eqs. (|59|) and l|62|l both constrain M to be less than approximately 
an inverse millimeter, or 10 -3 eV. It is remarkable that this is also the mass scale associated with the cosmological 



constant or dark energy causing the present acceleration of the universe 24 1. Granted, this is tiny compared to typical 
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particle physics scales, such as the weak or the Planck scale, and thus our potential (|56|l suffers from fine-tuning. 
Nevertheless, it is our hope that whatever mechanism suppresses the scale of the cosmological constant from its natural 
value of 10 18 GeV down to 10 -3 eV might also naturally account for the energy scale characterizing our potential. 

The above constraints can be translated into bounds on the range of chameleon-mediated interactions in the 
atmosphere (m~l m ), in the solar system (m^ 1 ) and on cosmological scales today (m^" 1 ). From Eq. (EH , these are 
given by 

m Pl 



m Pl 



»T, = V^M + -£rPoe^ o/Mp ' , (63) 
1u pi 



where po w 10 29 g/cm 3 is the current energy density of the universe and <po is the corresponding value of 4> on 
cosmological scales. Substituting the above bounds on M, it is straightforward to show that, for n < 2 and (3 of order 
unity, 



m atm < 1 mm - 1 cm 



mg 1 < 10 - 10 4 AU 



mp 1 < 0.1 - 10 3 pc, (64) 

where the numbers on the right-hand depend somewhat on the value of n and /? (recall that 1 AU ~ 1.5 • 10 13 cm and 
lpc«3-10 18 cm). 

Thus, while (^-interactions are short-range in the atmosphere (short compared to the size of the atmosphere), they 
are rather long-range in the solar system. In particular, it is possible for the scalar field to be essentially free on solar 
system scales. This striking difference in behavior of the field in space compared to Earth is an original ingredient of 
our scenario and, as we will show in Sec. lIXI can lead to unexpectedly large signals for EP and fifth force experiments 
to be performed in orbit in the near future. On cosmological scales, we see that a power-law potential implies an 
interaction range which is smaller than the present size of the observable universe, H^ 1 ~ 10 9 pc. It follows that mo 
is too large for 4> to be rolling on cosmological time scales today. A general study of potentials and their viability for 
quintessence models of the universe will be presented in a future paper 



VII. SOLAR SYSTEM TESTS OF GRAVITY 



In this Section, we discuss the constraints on scalar-tensor theories from planetary orbits, in particular from solar 
system tests of the EP and fifth force fSec. lVIIA^ . as well as post- Newtonian corrections fSec. lVH"*B|l . In short, these 
constraints are all satisfied because large bodies, such as the Sun and the Earth, all have a thin shell which greatly 
suppresses the 0-force between them. 
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A. Solar System Tests of EP and Fifth Force 

Precise measurements of the lunar orbit from laser ranging constrains the difference in free-fall acceleration of the 
Moon and the Earth towards the Sun to be less than approximately one part in 10 13 |IJ. That is, denoting their 
respective acceleration by aMoon and a® , we have 

\a M oon-a e \ .£ 10 -i 3) (65) 



a N 

where a at is the Newtonian acceleration. 

In our model, this difference is naturally very small since the Sun, Earth and Moon are all subject to the thin shell 
effect. We have already imposed that the Earth (and its atmosphere) have a thin shell. So must the Sun, therefore, 
since its Newtonian potential is larger than that of the Earth. Hence we only need to show that the same holds true 
for the Moon. But, assuming 4>c ^> (pMoom this trivially follows from Eq. 

A^Moon ^ < 1CT 5 , (66) 

^ Moon 

where we have used $® = 10~ 9 and &Moon = 1CP 11 . 

Hence the 4> profile outside each of these bodies is given by Eq. I|27[) with ~ m,Q and </>oo = <pc- Assuming 
m G } > 1 AU, since this yields maximal violation of the EP, it is then straightforward to show that the acceleration 
of the Earth towards the Sun is given by 



- = " ■ ■ < 1 + 18,3 l^f J {-W ) * "~ ■ 1 + m 1=5? J t • (67) 



32 /Ai?©\ fAR Q \\ J 2 f AR e \ 2 $ e 



while for the Moon 



aMoon = a N ■ < 1 + 18/3 — — — } w a N ■ < 1 + 18/3 — — - w } . (68) 



Substituting <j> = 10 , = 10 and $Moon = 10 , this gives a difference in free-fall acceleration of 

\a M oon -a®\ Kf3 2( A^eV < ^ , 1Q -i4 (6Q) 
a N \ i?® J 

where we have used Eq. (|34|l in the last step. This satisfies the bound in Eq. 1)65(1 for reasonable values of (3. 

We next consider solar-system tests of the existence of a fifth force. Deviations from a l/r 2 force law, for instance 
due to the exponential factor in Eq. I|44|l . contribute an anomalous component to the perihelion precession of planetary 
orbits in comparison with the predictions of GR. For instance, Lunar laser-ranging measurements lead to the constraint 
a < 10~ 10 for a fifth-force with range A ~ 10 8 m 34]. A similar analysis for the orbits of Mercury and Mars gives 
a < 10~ 9 for the range A ~ 1 AU [3J|. (Not surprisingly, these tests are most sensitive to a fifth force whose range is of 
the order of the distance between the Sun and the orbiting body.) In our model, these celestial objects are all subject 
to the thin-shell effect, and, just as with the EP analysis above, the screening mechanism makes the constraints from 
perihelion precession trivial to satisfy. 



19 



B. Tests of Post-Newtonian Gravity 

To estimate the constraints from post- Newtonian corrections, consider the 0-profile due to the Earth given by 
Eq. (33: 

where we have neglected the exponential factor. Comparison with the expected profile if there were no thin-shell 
suppression, given by Eq. J20J), we see that the exterior solution above corresponds to that of a massless scalar with 
effective coupling 

A /? 

Pef f = 3[3— ^<3/3-10" 7 , (71) 
-ft© 

where in the last step we have used the condition that the atmosphere has a thin-shell (see Eq. (|34() '). Treating our 
model as a Brans-Dicke theory with effective coupling constant /3 e // given above, which is a good approximation 
in the solar system since the chameleon behaves essentially as a free field, it is straightforward to show that the 
corresponding effective Brans-Dicke parameter, u>bd, is given by Q 

3 + 2lu bd = — l — > 6 ■ 10 12 /T 2 . (72) 

The tightest constraint on Brans-Dicke theories comes from light-deflection measurements using very-long-baseline 
radio interferometer 8]: lobd > 3500. We see that this is easily satisfied in our model. Similarly, one can show that 
the constraint from the decay of the orbital period of binary pulsars, ujbd > 100, is trivially satisfied. (Note that 
density-dependent effective couplings were previously noted in a different context [33j.) 

VIII. TESTS OF THE STRONG EP 

Our discussion of the EP has so far been restricted to the weak EP, which essentially states that the laws of gravity 
are the same in any inertial frame. General Relativity, however, satisfies a stronger version of the EP in the sense 
that all laws of physics, including non-gravitational interactions, assume the same form in any inertial frame (local 
Lorentz invariance), and the various parameters describing these non-gravitational forces, such as the fine-structure 
constant, are independent of space and time (local position invariance). 

At any space-time event, one can find a coordinate system in which the Einstein- frame metric equals the 
Minkowski metric r\u,v- Since all other metrics g^l appearing in the action Q are conformally-related to g^ v , all g^l 
are proportional to rj^ in this frame. All laws of physics are therefore Lorentz invariant at that space-time event, 
and thus our model satisfies local Lorentz invariance. 

n 

The most stringent constraint on spatial variations of couplings comes from the Vessot-Levine experiment _23] which 
measured the redshift between a hydrogen-maser clock flown at an altitude of 10 4 km and another one on the ground. 
In a theory which has LPI, as in GR, the redshift z is given by the difference in Newtonian potential A$ between the 
emitter and the receiver j^l- As shown in Will LPI violations generate an extra contribution Az to the redshift, 
of the form 



Az = 7 ■ A<f> , 



(73) 
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where 7 is a constant that depends on the Newtonian potential of the emitter, with 7 = corresponding to the case 
of no LPI violation. The bound from the Vessot-Levine experiment is I7I < 10 -4 . 

To estimate 7 in our case, recall from Eq. @ that a test particle of the matter field ipm follows geodesies of the 
metric gjfl related to the Einstein- frame metric g^ by g$ — e 2 ^/ Mpi g^. Thus, a constant mass scale mfi in 
the ip m -frame is related to a (^-dependent mass scale m(<f>) in Einstein frame by the rescaling m((f>) — e^^ Mp 'm^\ 
Similarly, a i/^-clock with frequency i/W is measured in Einstein frame to have a frequency v(<fr) — e^ i ^/ Mpl v^ l \ 

Assuming j3(p/Mpi <C 1, the ^-dependence of v therefore yields the following extra contribution to the redshift: 

Az^--^-((fr(r em )-q)(r rec )), (74) 
Mpi 

where we have dropped the superscript (i), and where r ern (r rec ) is the distance between the emitter (receiver) and 
the center of the Earth. For the Vessot-Levine experiment, we have r em « 10 km « 2_R^ and r rec > Rq, and it 
follows from Eqs. I|36|) that <p{r em ) rs <Pq and 4>{r rec ) s» <p a tm, respectively. Since A$ s» $e/2 in this case, massaging 
Eq. gives 



Az = -I20 z [ I" A$, (75) 

from which we can read off the corresponding value of 7: 

| 7 |= 12/32^© </j2. 10 -6_ (76) 

This comfortably satisfies the Vessot-Levine bound of I7I < 10~ 4 . 

Time-variation of coupling constants are constrained by geophysical measurements (such as the Oklo Nuclear 
Reactor 1^3), by the study of absorption lines in quasar spectra and by nucleosynthesis. Recent analysis indicates 
that the fine-structure constant has evolved by more than one part in 10 5 over the cosmological redshift range 
0.2 < z < 3.7 |5j], thus suggesting that LPI does not hold in the Universe. 

In our scenario, the various coupling constants are determined by <f>, whose value depends on the local density. 
Thus, even though the coupling constants may vary on cosmological scales today, the fact that the density of the 
Earth is constant in time implies that the various couplings measured on Earth do not vary significantly. In particular, 
this implies that constraints from the Oklo Nuclear Reactor are easily evaded in our model. (See [2S| for a related 
discussion of a density-dependent fine-structure constant.) 

In any case, we find that the time-variation of coupling constants and masses on cosmological scales is very small 
in the case of the power-law potential of Eq. 1)56(1 . To see this, consider once again a constant mass scale in 
the matter frame, with corresponding m((f>) — e^ i ^/ Mpl m^ in the Einstein frame. Hence, the time variation of m((j>) 
between nucleosynthesis and the present epoch, say, is given by 

Ar 



'' TT - ^ _ ^bbn) , (77) 
Mpi 

where 4>bbn is the value of <p at nucleosynthesis. For the power-law potential considered earlier, V{<p) — M 4+ ™0 _ ", 
the bound of M < (1 mm) -1 derived in Sec. I VII gives 

— < /3 • 10- 11 , (78) 
m 

which satisfies all current astrophysical and cosmological bounds. In particular, it appears that the power-law potential 
cannot account for the variation of the fine-structure constant reported by Webb et al. Q ■ A more detailed analysis of 
time- variation of coupling constants in our model, including more general potentials, will be presented elsewhere [2^. 
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IX. NEW PREDICTIONS FOR NEAR-FUTURE SATELLITE EXPERIMENTS 

Our scenario has the remarkable feature that the physical characteristics of the scalar field can be very different in 
the laboratory than in space. We have seen, for instance, that the range of the interactions it mediates is of order 
1 mm in the atmosphere while being greater than 10 AU in the solar system. Hence, as we will show in this Section, 
the predicted strength of EP violations and fifth force can be strikingly different in orbit than in the laboratory. This 
is particularly interesting as this decade should witness the launch of several satellite experiments that aim at making 
these measurements. Three of them, namely STEP, GG and MICROSCOPE, will test the universality of free-fall 
with respective expected sensitivity of 10~ 18 , 10~ 17 and 10 -15 for the Eotvos parameter rj. Another satellite, the SEE 
Project, will further attempt to measure a fifth force between two test bodies in orbit. 

In this Section we will show that, for a wide range of parameters, our model predicts that these satellite experiments 
will see a strong signal. It is in fact likely that the signal will be stronger than previously thought possible based on 
laboratory measurements. For instance, the SEE capsule could detect corrections to the effective Newton's constant 
(including fifth force contribution) of order unity compared to the value measured on Earth or inferred from planetary 
motion. Moreover, STEP, GG and MICROSCOPE could measure violations of the EP with r] > 10~ 13 , which is larger 
than the current bound from the (ground-based) Eot-Wash experiment. Such outcomes would constitute a smoking 
gun for our model, for it would otherwise be difficult to reconcile the results in space with those on Earth. 

Let us begin with the SEE satellite, which will orbit the Earth at an altitude of approximately 10 3 km. This multi- 
faceted experiment will test for deviations from 1/r 2 in the force law, search for violations of the EP and measure 
Newton's constant G to one part in 10 7 . This will be achieved by accurately determining the orbit of two tests masses 
as they interact gravitationally with each other and with the Earth. For our purposes we shall focus on the expected 
value of G, including the contribution from the fifth force mediated by <f>. 

We first show that it is possible for the satellite not to have a thin shell. This is desirable in this case in order to 
maximize the signal. That is, we derive the conditions under which ARsee/Rsee > 1 holds for the SEE capsule. 
According to the current design, the satellite will have a total mass of ~ 2000 kg and an effective radius of Rsee ~2m 
(although it has cylindrical symmetry, we approximate the capsule as a sphere of equal volume for simplicity). Thus 
its Newtonian potential is &see ~ 10~ 24 sw 10~ 15 $&. Moreover, at an altitude of 1000 km, the chameleon field 
assumes the value 4>(Tsee) ~ 4>G, as seen from Eq. JSSJ. Therefore, the condition ARsee/Rsee > 1 requires 

^® > 10- 15 . (79) 

Combining this with the condition that the atmosphere have a thin shell, Eq. (|57|) , we find that the allowed range for 
which the SEE satellite will not have a thin shell is 

lO-is < ^® < 10- 7 . (80) 
-fie 

This constitutes a wide region in parameter space. Following the analysis of Sec. lVIl for a potential V(<j>) = M 4+n (f)~ n 
with n = 1/3 and [3 = 10, this requires the range of interactions in the atmosphere, m~ tm , to be anywhere between 
0.1 /itm and 1 mm. 

Therefore, if the inequality in Eq. I|79|l is satisfied, the scalar field 4> is only slightly perturbed by the satellite, and 
thus the field value within the capsule is <fi ~ <pa- Moreover, since tUq 1 is much larger than the size of the satellite, <fi 
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behaves as a massless field. Hence, if the test bodies have mass Mj and coupling i = 1, 2, they will experience a 
total force, gravitational plus ^-mediated, given by 

\F\ = (1 + 2ftA) , (81) 

where we have restored Newton's constant G = (87rAf|, ; ) _1 . This implies, therefore, that the SEE satellite will 
measure an effective Newton's constant, G e // = G(l + 2/3i/32), which differs by order unity from Newton's constant 
measured on Earth. 

Since this is an important prediction, let us summarize how it was obtained. First of all, the key ingredient is that 
the SEE experiment is performed in orbit, where the range of ^-mediated interactions is greater than 10 AU. Thus, 
small perturbations in <f> are essentially massless on the scale of the capsule, and therefore the force law between 
two test bodies is proportional to 1/r 2 to a very good approximation. A second crucial condition is that there is no 
thin-shell effect within the satellite, as expressed in Eq. Q79JI . This ensures that there is no suppression of the fifth 
force, and the coupling strength is therefore of order 1 . Hence, for (3 of order unity, the effective Newton's constant 
receives large corrections from the fifth force mediated by <f>. As mentioned earlier, such an outcome would constitute 
strong evidence that a mechanism of the kind we are proposing is at play in Nature. 

Our model also has surprising implications for satellite experiments that aim at measuring EP violations in orbit, 
such as MICROSCOPE, GG and STEP. These experiments will attempt to measure a difference in free-fall acceleration 
between two concentric cylinders of different composition with expected sensitivity of one part in 10 15 , 10 17 and 10 18 , 
respectively. 

We focus on the STEP satellite for concreteness. Although Eq. (|79|l pertains to the SEE capsule, a similar condition 
is obtained for STEP, since its physical characteristics are not too different from those of SEE. Hence, there is no 
thin-shell effect if Eq. (j79(l holds, and the 4> profile within the satellite is well-approximated by Eq. I|35|l . It is then 
straightforward to show that the Eotvos parameter for Be and Nb test cylinders, as appropriate for STEP, is given 
by 

r, « 10- 4 /3 2 ^. (82) 

Combining with Eq. (|80|l . we find the allowed range 

f3 2 ■ 10- 19 < r\ < (3 2 • KT 11 , (83) 

which falls almost entirely within the range of sensitivity of STEP. Moreover, we see that (3 can be larger than 10~ 13 , 
the current bound from experiments performed on Earth [9j. In other words, it is possible that MICROSCOPE, GG 
and STEP will measure violations of EP stronger than currently thought to be allowed by laboratory measurements. 
This is another striking manifestation of the fact that scalar field dynamics in our model are very different on Earth 
than in orbit. 



X. DISCUSSION 



In this paper, we have presented a novel scenario in which both the strength and the range of interactions mediated 
by a scalar field <f> depend sensitively on the surrounding environment. In a region of high density, such as on Earth, 
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the mass of the scalar field is sufficiently large, typically of order 1 mm -1 , to evade constraints on EP violation from 
laboratory experiments. Meanwhile, the field can be essentially free on solar system scales, with a typical Compton 
wavelength of 100 AU. 

We have argued that our scenario satisfies all existing bounds from ground-based and solar-system tests of gravity. 
It will be left for future work to study our model in the context of strongly gravitating systems, such as black holes [40j. 
Implications for cosmology will be analyzed in a separate paper 24] , where we will show in detail that chameleons are 
consistent with cosmological constraints on the existence of non-minimally coupled scalars, such as the bound on the 
time-variation of G from nucleosynthesis, for example. Moreover, we will describe how the field dynamically reaches 
the minimum of the effective potential, an element that was assumed a priori in the present work, and show how our 
scenario naturally provides a solution to the old moduli problem [-111 ]. 

The striking difference in interaction range in the solar system versus on Earth can lead to unexpected outcomes 
for experiments that test gravity in space. It was argued that the SEE Project could measure an effective Newton's 
constant drastically different than that on Earth. Meanwhile, the MICROSCOPE, GG and STEP satellites could 
detect violations of the EP stronger than currently allowed by ground-based experiments. Either outcome would 
effectively constitute a proof of the existence of chameleons. More importantly, these possible surprises for space 
experiments will hopefully strengthen the case for launching these satellites. 

Also, the door is open for conceiving of new table-top experiments to probe for the existence of chameleons, for 
instance by exploiting the fact that the properties of chameleons are sensitive to the surrounding matter density. For 
example, one could imagine doing atomic physics experiments inside or in the vicinity of a massive oscillating shell of 
matter. These oscillations would induce a time-dependence in the profile of the chameleon and thus in the frequency 
of emission lines 0|. It would also be interesting to investigate the behavior of the field at intermediate altitude, 
for instance at about 40 km where balloon experiments take place. This would require a more realistic modeling of 
the atmosphere than presented here. Such an analysis could reveal that it may be possible to detect chameleons by 
performing gravity experiments in balloons . 
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